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Abstract 

The Beam Energy Scan (BES) program is being pursued at RHIC to study the QCD phase di- 
agram, and search for the possible QCD phase boundary and possible QCD critical point. The 
data for Phase-I of the BES program have been collected for Au+Au collisions at center-of-mass 
energies ( yfsm) of 7.7, 1 1.5, 19.6, 27, and 39 GeV. These collision energies allowed the STAR 
experiment to cover a wide range of baryon chemical potential fig (100-400 MeV) in the QCD 
phase diagram. We report on several interesting results from the BES Phase-I covering the high 
net-baryon density region. These results shed light on particle production mechanism and freeze- 
out conditions, first-order phase transition and "turn-off" of QGP signatures, and existence of a 
critical point in the phase diagram. Finally, we give an outlook for the future BES Phase-II 
program and a possible fixed target program at STAR. 



1. Introduction 

The results from top RHIC energies suggest the existence of the Quark Gluon Plasma (QGP) {3] . 
The main task at hand now is to study the properties of the QGP and establish the QCD phase 
diagram. Lattice QCD calculations predict the transition between QGP and the hadronic gas as 
crossover at /j.b = 0, while at large /jb they predict a first order phase transition yfl. A point where 
the first order phase transition ends is called the critical point. Experimentally, the QCD phase 
diagram can be studied by colliding heavy ions at varying beam energies that can provide a T-/j.b 
region for each energy. Then, one can look at the various signatures for the phase boundary and 
critical point. We present results on (a) freeze-out parameters to get insight into the QCD phase 
diagram, (b) first and second coefficients (vi, V2) of the Fourier expansion of the angular dis- 
tributions, dynamical charge correlations, and nuclear modification factor Rep for the search of 
first-order phase transition and "turn-off" of QGP signatures, and (c) fluctuation measurements 
for the search of QCD critical point. 

The first STAR proposal for the BES program jj] was made in the year 2008. This was 
followed by a successful data taking and physics analysis of a test run below injection energies 
at ysAw = 9.2 GeV Jj]. The first phase of the BES program was started in the year 2010 with 
data taking at three low energies of 7.7, 11.5, and 39 GeV. In the year 2011, two more energy 
points were added at yfsNN = 19.6 and 27 GeV. 
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The results presented here are based on the data collected by the STAR detector. STAR 
covers a large acceptance of 2n in azimuth (0) and -1 to 1 in pseudorapidity (rj). One of the im- 
portant advantages which STAR possesses for the BES program is its almost uniform acceptance 
for different identified particles and collision energies at midrapidity. The main tracking device 
at STAR is the Time Projection Chamber (TPC), which provides momentum as well as particle 
identification (PID). For higher transverse momentum (pr) region, the Time Of Flight (TOF) 
detector is quite effective in distinguishing different particles. Particles are identified using the 
ionization energy loss in TPC and time-of-flight information from TOF [5]. The centrality se- 
lection is done in STAR using the uncorrected charge particle multiplicity measured in the TPC 
within I77I < 0.5 H. 

2. Accessing QCD Phase Diagram 

The QCD phase diagram is the variation of temperature T and baryon chemical potential ps- 
These quantities can be extracted from the measured hadron yields. Transverse momentum spec- 
tra for the BES Phase-I energies are obtained for n, K, p, A, E, K^, and <p (g]. The particle ratios 
are used to obtain the chemical freeze-out (a state when the yields of particles get fixed) condi- 
tions using the statistical thermal model (THERMUS) [7]. The two main extracted parameters 
are chemical freeze-out temperature T c h and fig- Figure Q] (left panel) shows the variation of the 
extracted chemical freeze-out parameters using the Grand-Canonical Ensemble (GCE) approach 
of THERMUS for different energies and centralities [8]. We observe that at lower energies, T c h 
shows a variation with /ig as a function of centrality. Both T c h and ps values seem to be de- 
creasing from central to peripheral collisions. It may be noted that use of Strangeness Canonical 
Ensemble (SCE) shows an opposite trend, i.e. the T c h values seem to increase from central to 
peripheral collisions |8p. However, it is found that the x 2 /NDF in case of SCE is higher for 
peripheral collisions [8]. More investigations are ongoing. In GCE the energy and quantum 
numbers or particle numbers are conserved on average through the temperature and chemical 
potentials. GCE is widely used in heavy-ion collisions. In the SCE, the strangeness (S) in the 
system is fixed exactly by its initial value of S, while the baryon and charge contents are treated 
grand-canonically. Although the T c h shows an opposite trend between GCE and SCE, both re- 
sults suggest the variation of chemical freeze-out parameters with centrality at lower energies. 

The particle spectra can be used to obtain the kinetic freeze-out (a state when the spectral 
shapes of particles get fixed) conditions using the Blast Wave (BW) model [11]. The BW model 
is used to simultaneously fit the n, K, p spectra and the two relevant extracted parameters are 
kinetic freeze-out temperature T^m and average flow velocity (J3). Figure[T](right panel) shows the 
variation of kinetic freeze-out parameters for different energies and centralities |8J]. We observe 
that at a given collision energy, there is an anti-correlation between 7\j n and (J3). For a given 
collision centrality, the freeze-out temperature at high energy is lower and the average collectivity 
velocity (fi) is larger due to expansion. 

3. Search for First Order Phase Transition & Turn-off of QGP Signatures 

3.1. Directed Flow 

The directed flow vi is calculated as (cos(0 - where and W\ are the azimuthal angle 
of the produced particles and orientation of the first-order event plane, respectively. The directed 
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Figure 1: (Color online) Left panel: Variation of T c t with un for different energies and centrali- 
ties. The curves represent the theoretical calculations j^ fioll . Right panel: Variation of 7\i n with 
(fi) for different energies and centralities. Errors in both panels represent the quadrature sum of 
systematic and statistical errors. 



flow measurements near midrapidity for protons are considered as sensitive to the equation of 
state (EOS) and hence can be considered a phase transition signal [12]. Figure [2] shows the vi 
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Figure 2: (Color online) Directed flow slope 
(dvi/dy', /=y/ybeam) for n ', p, p, and net- 
protons {p-p) near midrapidity as a func- 
tion of beam energy for mid-central (10-40%) 
Au+Au collisions. The shaded band refers to 
the systematic uncertainty on net-proton mea- 
surements. 




Figure 3: (Color online) The difference in 
v'2 between particles and their corresponding 
anti-particles as a function of beam energy in 
0-80% Au+Au collisions. The curves repre- 
sent fits to data points as discussed in text. 
Both statistical (vertical lines) and systematic 
errors (caps) are shown. 



slope (dvi/dy', where y' - y/ybeam and y is rapidity), near midrapidity as a function of beam 

i — [ 

energy for the mid-central (10^-0%) Au+Au collisions [13]. The pion and anti-proton v\ slopes 
show negative values for all the beam energies studied. The proton v\ slope changes sign while 
going from 7.7 GeV to 11.5 GeV and then stays negative up to 200 GeV. However, the net- 
protons vi slope (obtained using Vj slopes of p, p and ratio of p/ p) changes sign two times going 
from lower to higher energy while showing a dip around yfsNN- 10-20 GeV. More studies are 
needed in order to understand these interesting observations. 
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3.2. Elliptic Flow 

The elliptic flow v% is calculated as (cos 2(0 - ^2)), where T^ is orientation of the second- 
order event plane. Elliptic flow mainly probes the early stages of heavy-ion collisions. At top 
RHIC energy of 200 GeV in Au+Au collisions, the elliptic flow scaled by the number of con- 
stituent quarks (n q ) vs. {mj - mo)/n q (where mj - Jp^, + mfy shows a scaling behavior where 
mesons and baryons have similar values at intermediate pj. This is referred to as the number of 
constituent quark (NCQ) scaling JI4I1 . It is an established signature of partonic matter formed in 
Au+Au collisions at 200 GeV and deviations from such scaling would indicate the dominance 
of hadronic interactions. Hence breaking of NCQ scaling at lower energies could be an indi- 
cation of a "turn-off" of QGP signatures. Figure [3] shows the difference in V2 of particles and 
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Figure 4: (Color online) Vz/n q as a function of (raj - mo)/n q for different particles in Au+Au 
collisions at yfsim- 7.7, 1 1.5, 19.6, 27, 39 and 62.4 GeV. The errors shown are statistical only. 

corresponding anti-particles as a function of beam energy The curves represent fits to data 
points with functional form: /a V2 ( s/snn) = a tJsnn +c. The V2 difference between particles and 
anti-particles is observed to increase when we go towards the lower energies. At low energies, 
v 2(n ) > V2(tt + ), V2{K + ) > V2(K~), and V2(baryons) > V2(anti-baryons). This difference between 
particles and anti-particles suggests that the NCQ scaling among particles and anti-particles is 
broken. However, the observed difference between V2 of particles and anti-particles could be 
qualitatively explained by the models incorporating baryon transport at midrapidity and hadronic 
interactions [ 16]. We also observe that the baryons-mesons splitting for V2 vs. mj - mo starts to 
disappear for anti-particles at 1 1 .5 GeV and below. Figure|4]shows the \'oJn q vs. (mj-mo)ln q for 



different particles for 7.7-62.4 GeV B15I1 . We observe that results for all the particles are 

consistent among each other within ±10% level, except for the 0-mesons. At the largest mj - mo 
the 0-meson data points deviate by 1.8<x and 2.3cr for ysiviv = 7.7 and 11.5 GeV, respectively. 
Since ^-mesons have smaller hadronic interaction cross-section, their smaller v'2 values could in- 
dicate that the hadronic interactions start to dominate over partonic effects for the systems formed 
at beam energies below y^wv = 

11.5 GeV HI. However, as can been seen from the figure, a 
higher statistics data are needed to extend the mj - mo range and significance of the deviation 
observed. 
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3.3. Dynamical Charge Correlations 

The dynamical charge correlations are studied through a three-particle mixed harmonics az- 
imuthal correlator 11811 . y - (cos(<f> a + 4>p - 2 v Frp)). This observable represents the difference 
between azimuthal correlations projected onto the direction of the angular momentum vector and 
correlations projected onto the collision reaction plane. It is suggested that the difference in the 
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Figure 5: (Color online) (cos((f> a + (f>p — 2 v Prp)} as a function of centrality for Au+Au collisions at 
y*AW =7.7-200 GeV. For comparison, results for Pb+Pb collisions at 2.76 TeV are also shown. 
Errors are statistical only. 



correlations between same sign and opposite sign charges in heavy-ion collisions could be related 
to local parity violation if there is a deconfinement and a chiral phase transition lfl9ll . This is also 
referred to as Chiral Magnetic Effect (CME). At top RHIC energies, we observed a separation 
between the correlations of same and opposite sign charges. If this difference can be attributed to 
the QCD phase transitions, the absence of such observation could be an indication of the system 
which did not undergo the phase transition. Hence, the observable could be useful to locate the 
energy in the BES program where the QGP signature "turns off". Figure [5] shows the results for 
the beam energies from 7.7-200 GeV as a function of centrality 12011 . For comparison, Pb-Pb 
results from ALICE are also shown 112 ill which are observed to be consistent with the results 
from top RHIC energy. The separation between same and opposite sign charges decreases with 
decreasing energy and vanishes below y/s^N =1 1.5 GeV. 

3.4. Nuclear Modification Factor 

Another established signature of QGP at top RHIC energy is the Nuclear Modification Factor 
Rcp, which is defined as ratio of yields at central collisions to peripheral collisions, scaled by the 
corresponding number of binary collisions. It has been observed that at high pj, the Rqp of 
different particles is less than unity [22], which is attributed to the energy loss of the partons in 
the dense medium. In the absence of dense medium, there may not be suppression of high pj 
particles, which can serve as an indication of "turn-off" of a QGP signature. Figure [6] shows the 
Rcp of various strange hadrons such as K®, <p, A, H~, and Q. in Au+Au collisions at -\/Jnn =7.7- 
39 GeV Q. We observe that for p T > 2 GeV/c, the R C p(K°) is less than unity at 39 GeV and 
then the value increases as the beam energy decreases. For ysAW < 19.6 GeV, Rcp(K^) is 
above unity, indicating decreasing partonic effects at lower energies. For (un)-identified charged 
hadrons results, please refer to [24]. 

5 



'"I""!""!""!'" 

39 GeV 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

27 GeV 

STAR Preliminary 




in i |i ■ ■ ■ | n n | ■ ■ ■ i| n n 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

11.5 GeV 



i r i 1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

7.7 GeV 



+4i 



01234 1234 1234 1234 12345 

P T (GeV/c) 

Figure 6: (Color online) R C p (0-5% / 40%-60%) for A, H~, Q and R CP (0-10% / 40%- 
60%) for 0-meson in Au+Au collisions at y/l^N -1.1-39 GeV. Errors are statistical only. Grey 
bands represent the normalization error from AW- 



4. Search for QCD Critical Point 

If a system passes close to a critical point, the correlation length £ is expected to diverge. 
The higher moments (skewness S and kurtosis k) of distributions of conserved quantities such 
as net-baryons, net-charge, and net-strangeness, have a better sensitivity to the correlation length 
compared to variance cr 2 B25I1 . Also, the moment products such as kct 2 and Scr can be related to 
the ratios of order susceptibilities calculated in Lattice QCD and HRG model as kct 2 = X'b IXb* 
and Scr — Xb* Ix^r lES- One of the advantages of using these products is that they cancel the 
volume effects. Figure [7] shows the kct 2 and S cr for net-protons as a function of beam energy for 
different collision centralities B27I1 . For comparison, the results are shown for Poisson expecta- 



tions and UrQMD model calculations [28]. The bottom panel shows the Scr values normalized 
by the corresponding Poisson expectations. We observe that the moment products show similar 
values for central 0-5% and peripheral collisions for yfsJiN =39-200 GeV. For beam energies 
below 39 GeV, they have different values for central and peripheral collisions. Their values are 
below Poisson expectations for ys/wv > 7.7 GeV for 0-5% central collisions. However, for 
peripheral collisions, moment products values are greater than Poisson expectations below 19.6 
GeV. The UrQMD model calculations show a smooth monotonic behavior as a function of col- 
lision energy. It may be noted that we have large uncertainties for data points below 19.6 GeV 
that calls for higher statistics data for these energies. For net-charge and particle ratio fluctuation 
results, please refer to 12^1 . 



5. The BES Phase-II and STAR in Fixed Target Mode 

The BES Phase-I program from STAR has several interesting results as discussed above. We 
have collected sizable data for Au+Au collisions over a wide range of collisions energies from 
7.7 to 39 GeV. However, to consolidate the findings from Phase-I, we need higher statistics at 
lower energies, especially at 7.7 and 1 1.5 GeV. As we have discussed above, statistics for several 
important observables, such as 0-meson V2 and higher moments of net-protons distributions, are 
not sufficient to draw quantitative conclusions. Also, there is more than a 100 MeV gap in fig 
between 11.5 and 19.6 GeV, so having one more energy point around 15 GeV would allow to 
confirm the trends happening between 1 1.5 and 19.6 GeV. 
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Table 1 : The collider mode y/s^N, fixed target 
mode -\/snn, and corresponding fig values. The 
fis values are listed for the central collisions cor- 
responding to fixed target yfsim 
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Figure 7: (Color online) ko 2 , Scr and Scr val- 
ues normalized by the Poisson expectation, as 
a function of collision energy. Vertical bars are 
statistical and caps are systematic errors. 



In view of these, STAR has proposed a BES Phase-II program. For this STAR has requested 
an electron cooling device at RHIC to increase the luminosity for collisions below yfsim = 20 
GeV. Simulation results indicate that with electron cooling, the luminosity could be increased by 
a factor of about 3-5 at 7.7 GeV and about 10 around 20 GeV [30]. An additional improvement 
in luminosity may be possible by operating with longer bunches at the space-charge limit in the 
collider l5Hl . So overall improvement in luminosity could be about 10-fold. The luminosity 
improvement will not only allow the precision measurements of the above important observables 
but will also be helpful in the measurements of rare probes such as dilepton production and 
hypertriton measurements [ 32] . To maximize the use of collisions provided at STAR for the BES 
program, we have proposed to run STAR in "fixed-target mode". The main motivation for this 
proposal is to extend the ^b coverage currently from 400 MeV to about 800 MeV so as to cover 
a large portion of the phase diagram. For this, we plan to install a fixed Au target in the beam 
pipe to perform the Au(beam)-Au(target) collisions. The data taking can be done concurrently 
during the normal RHIC running. This proposal will not affect the normal RHIC operations. 
Table[T]lists the proposed collision energies, corresponding fixed target center-of-mass energies, 
and baryon chemical potential values. 

6. Summary 

We have presented several interesting results from the STAR BES Phase-I program. Cur- 
rently, STAR covers a large range of {1b (20-400 MeV) in the phase diagram. We have observed 
a centrality dependence of freeze-out parameters at lower energies. We have looked at various 
observables for the search of first order phase transition, "turn-off" of the QGP signatures, and 
search of a QCD critical point. The observables such as v\ slope, Re? and V2 of PID hadrons, and 

higher moments of net-proton distributions, show some dramatic changes for the energy range 

? 



ysjvjv < 20 GeV. We have observed that proton V\ slope changes sign between 7.7 and 1 1.5 GeV. 
The net-proton vi slope changes sign twice as a function of beam energy. We have observed a 
difference in V2 of particles and corresponding anti-particles which increases with decreasing col- 
lision energy. The 0-meson V2 deviates from that of other particles, the charge separation signal 
vanishes, and the Rqp of is greater than unity for beam energy of 1 1 .5 GeV and below. The 
moment products have values below the Poisson expectation for 0-5% central Au+Au collisions. 
For peripheral collisions, they are above unity at y/s^N < 19.6 GeV. An electron cooling device 
at RHIC will significantly increase luminosity at beam energies below 20 GeV allowing the pre- 
cision measurements for several important observables for phase structure studies. In addition, 
we propose to run STAR in fixed target mode which will extend the fig coverage from STAR up 
to 800 MeV in the phase diagram. 

We acknowledge the support from DOE for this research. 
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